function Bh_init=setPriorMean_RW(modelSpec,isRandomWalk)
%initialize prior on Minnesota-type prior with ones on main diagonal 

%unpack basics
n  =modelSpec.modelSize; 
nP =modelSpec.nBLPlags;


Bh_init=zeros(n*nP+1,n);

%populate diagonal
Bh_init(2:n+1,:)=diag(1.*isRandomWalk); %prior mean (coefficients)
end